The thermal width of heavy quarkonia moving in quark gluon plasma 
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The velocity dependence of the thermal width of heavy quarkonia traveling with respect to the 
quark gluon plasma is calculated up to the NLO in perturbative QCD. At the LO, the width 
decreases with increasing speed, whereas at the NLO it increases with a magnitude approximately 
proportional to the expectation value of the relative velocity between the quarkonium and a parton 
in thermal equilibrium. Such an asymptotic behavior is due to the NLO dissociation cross section 
converging to a nonvanishing value in the high energy limit. 
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I. INTRODUCTION 



Using arguments based on the contraction of the Debye 
screeninglength in quark gluon plasma (QGP), Matsui 
and Satz[l| suggested J/ip suppression to be a signature 
of the formation of QGP in the early stages of a heavy 
ion collision. Indeed measurements at past SPS data [2] 
showed nontrivial suppression patterns that could be con- 
sistent with the original prediction. However, recent lat- 
tice calculations show that J/ ip will survive past the crit- 
ical temperature T c for the phase transition 0, 0, S EI 0] 
up to about 1.6 T c , while Xc and ip 1 will dissolve above 
T c . These findings suggest that a possible mechanisms 
for J/ip suppression could be the disappearance of feed- 
backs from the \c and ip'i together with the hadronic 
matter effects, such as the the nuclear absorption, the in- 
teractions with comovers, and the shadowing effect 
Using a more precise determination of the cold nuclear 
matter absorption cross section of charmonium through 
p-A collisions [nj, it was found that such a feed down 
suppression scenario was indeed favored in results from 
semi-central Pb-Pb collisions fill and from central In-In 
collisions at 158 GeV/nucleon fl2j ]. 

Phenomenologically, the statistical hadronization 
model [H, ITil. flBl fig . [l7| applied to the charmonium 
production [18| appears to be compatible with RHIC 
data[20l|. In a kinetic model, the charmonium is produced 
in the whole temporal evolution of the QGP [2l|, (2^, [H| . 
Further tests at LHC will discriminate between various 
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pictures, and a unified picture is expected to emerge. 

But before a simplified picture of J / ip suppression can 
be adopted, detailed properties of J/ ip above T c have to 
be investigated. Unfortunately, the present lattice cal- 
culations based on Maximum Entropy Method still have 
poor resolution, and is not able to reliably determine the 
thermal width or possible mass shift above T c 0,01 . How- 
ever, there are several recent works that can supplement 
the lattice calculation. A recent QCD sum rule calcula- 
tion using the running coupling constant and the gluon 
condensates extracted from a recent lattice data shows 
that the width of J/jjj may be broadened or the mass 
reduced just over T c [241 ], which is consistent with results 
from AdS/QCD (25|. In another recent work, the Debye 
screening length in the J/ip moving with respect to QGP 
was calculated in AdS/QCD [26]. In addition, the spec- 
tral functions of heavy quarkonia were extracted from 
the imaginary part of Green functions obtained from po- 
tentials fitted to lattice data [13, [H, H^, [3(| ■ In another 
work, the thermal width of J/ip was investigated using 
perturbative QCD up to the NLO [H|. The LO pertur- 
bative QCD calculation for the heavy quarkonium disso- 
ciation was invented by Peskin and Bhanot [32, 33] more 
than 20 years ago. Later, one of us redrived the LO result 
using Bethe-Salpeter amplitude [34| ■ which was further 
used to calculate the NLO results [35|]. In that work, 
the bound state of quarkonium was described by Bethe- 
Salpeter amplitude, and the perturbative QCD method 
was applied to calculate the decay process. The binding 
energy of the quarkonium and its radius were obtained 
by solving the Schrodinger equation with the potential 
energy extracted from the lattice QCD calculation [36| . 
At the NLO, it was found that the thermal width of J ftp 
increases as temperature increases, while it decreases in 
the LO. Moreover, the total width was found to grow 
to more than 250 MeV at 1.4 T c , assuming the thermal 
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quark gluon masses to be around 400 MeV. If the thermal 
width is so large, the J/ifi initially formed at 1.6 T c , will 
be dissociated immediately and not be able to escape the 
quark gluon plasma until it cools down further to near 
T c . 

Another important aspect to be considered in a realis- 
tic heavy ion collision is the velocity of the J/ tp with re- 
spect to the QGP. In particular, at LHC, more energetic 
heavy quarkonia are expected to be produced and ab- 
sorbed. Previously the Debye screening length between 
a heavy quark and a heavy anti-quark pair moving with 
respect to the QGP was calculated in a kinetic theory 
approach 37f . It was found that the Debye screening 
length becomes shorter as the velocity increases, because 
the parton density enhances in the heavy quarkonium 
rest frame [371 \3M- Recently, the screening length was 
investigated in a AdS/CFT calculation (2(| [39J, where it 
was found to be approximately proportional to [1— v 2 } 1 ^ 4 , 
v being the velocity of J/ip with respect to the QGP. 
But whether AdS/CFT calculations represents real QCD 
phenomena still remains controversial. Therefore, in this 
work, we will extend a previous NLO perturbative QCD 
calculation for the thermal width of a quarkonium at rest 
[3l| to that of a quarkonium moving with respect to the 
QGP. The works mentioned above [26|, [13, [H, Hi| antici- 
pated the shortening of the Debye screening length when 
the heavy quark and anti-quark pair moves in the quark 
gluon matter or equivalently when the matter is boosted, 
which may be called as the change of static properties. 
On the other hand, our results provide another compo- 
nent of the change through the thermal width of the 
quarkonia as it moves through the matter. This may be 
called as the change of dynamical properties. The short- 
ening of Debye screening length means the radius and 
the binding energy of heavy quarkonia should change as it 
moves through the matter. However, in this work, the ra- 
dius and the binding energy of static quarkonia are used 
for the purpose of investigating the change of a purely 
dynamical property. We find that at the LO, the width 
decreases with increasing v, which is caused by the van- 
ishing of the dissociation cross section of the quarkonium 
by an energetic gluon. However, at the NLO the thermal 
width increases with v. This is due to the nonvanishing 



asymptotic cross section between the quarkonium and 
an energetic parton at the NLO, which for the coulomb 
wave function of the quarkonium scales as the square of 
the Bhor radius. In section II, we briefly review formu- 
las used throughout this paper. In section III, we apply 
these formulas to J/ip and T. Some discussions are given 
in section IV. In the Appendix, we derive the asymptotic 
form of the dissociation cross section of a quarkonium in 
the high energy limit. 

II. THERMAL WIDTH OF A HEAVY 
QUARKONIUM IN PQCD 

The width of a hadron in the vacuum comes from its 
spontaneous decay. On the other hand, its thermal width 
results from its interactions with the surrounding thermal 
particles. This thermal width of a quarkonium moving 
with velocity (3 in the medium is defined as 



d 3 k 



n(k )v re i((3)a(k,/3), 



(1) 



where d p is the degeneracy factor, n(k ) is the distri- 
bution function of the thermal particle, v re i the relative 
velocity between decaying particle and the thermal par- 
ticle, and cr(fc, ft) their energy-dependent elementary dis- 
sociation cross section. Because we are interested in the 
decays of heavy quarkonia in the quark-gluon plasma, the 
decaying particle is a heavy quarkonium and the thermal 
particles are light quarks and gluons. In this work, we 
considered only three light flavors. The dissociation cross 
section can be written as 



a(k,P) 



1 



Av rel (p)E^(p)E p (k) 



d(p.s.)\M\ 



(2) 



The first factor on the right side is the inverse of initial 
flux. E p (k) are the energies of a quarkonium and 

a parton respectively, p.s. means phase space of final 
states, and M is the spin-averaged invariant amplitude. 
The invariant amplitudes for the decay of a quarkonium 
by partons are listed in [U [35[. Substituting Eq. @ 
into Eq. (TTJ we find, 
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Here the positive sign in the denominator is for a 
fermion and the negative sign for a boson medium. 
The function 6 + means only positive energy is allowed, 
and u is the four velocity of the thermal bath, here- 
after called the Lab. frame. In the Lab. fame, u = 
(1,0,0,0), whereas in the quarkonium rest frame, u = 



(1 - /3 2 )-V2(i 5 0,0,-/3). In the second equality of Eq. 
([3]), all factors except E$(/3) are Lorentz invariant. In 
fact, because i?$ = 7m$, where 7 = 1/ yl — (3 2 and 



m$ is the mass of a quarkonium, r yTp 
invariant quantity, and we have 
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FIG. 1: The variation of the thermal widths V of J/tp in the 
LO (upper) and in the NLO(lower) as a function of their veloc- 
ity in QGP for various temperatures, and assuming the ther- 
mal mass of a parton to be 400 MeV(left) or 600 MeV(right) 



= r e// 

quarkonium rest frame' 



(4) 



where r^{ rkonium rcsj; frame is calculated in the rest frame 
of the quarkonium with a moving medium. This relation 
just reflects time dilation. 



III. THE THERMAL WIDTH OF MOVING J/ip 
AND T IN QGP 

We now apply the preceding formulas to calculate the 
thermal width of the J/tp and the T. All calculation will 
be performed in the Lab frame, and T will be used to 
represent T^JL from now on. The degeneracy factor d p 
in Eq. Jl} is set to 16 for the gluons and to 36 for the 
quarks. For numerical purpose in these calculations, the 
T c value is taken to be 170 MeV. The figures in the upper 
panel of Fig. [T] show the variation of the thermal width 



T oi J/tp in the LO of perturbative QCD as a function 
of the velocity of the quarkonium in QGP. To the LO, 
the elementary dissociation process is a thermal gluon 
dissociating the quarkonium into a cc pair. Since the 
gluons acquire an effective thermal mass in the QGP, we 
introduce a constant thermal mass of either 400 or 600 
MeV[S(|, respectively to represent a value at the lower 
and upper limit. The graphs in the left figure are ob- 
tained with a thermal gluon mass of 400 MeV, and the 
right figures with 600 MeV. The dissociation cross sec- 
tion in the LO has the maximum value when the absorbed 
energy is slightly above the threshold, and then rapidly 
decreases as the gluon energy increases further [34], [35| • 
However, since the thermal mass of an absorbed gluon is 
larger than the energy that gives the maximum cross sec- 
tion, larger thermal mass will give smaller cross section. 
The binding energy of a J/tf) in QGP is obtained by solv- 
ing the Schrodinger equation with a potential extracted 
from the lattice data, from which we find that it varies 
from 36.4 MeV to few KeV as the temperature changes 
from 1.13 T c to 1.65 T c HIS!]. When the thermal mass 
of the absorbed gluon becomes larger, the dissociation 
begins at higher energy and the thermal width becomes 
smaller. This is why the thermal width with thermal 
gluon mass of 600 MeV is smaller than that with 400 
MeV. The graphs at the LO also show that the thermal 
width decreases as the velocity of J/tp in QGP increases. 
The reason for it is again simple to understand. As can 
be seen from Eq. ([!} , the thermal width is a convolution 
of elementary cross section , the relative velocity between 
the J/tp and the colliding parton, and the thermal dis- 
tribution function of the parton. However, as J/ip moves 
across the parton matter, many slow partons becomes 
fast partons as seen by the J/tp- As the dissociation cross 
section by fast parton is small in the LO, due to the small 
wave function overlap, the thermal width decreases. 

The figures in the lower panel of Fig. [1] show the ther- 
mal width of J/tp as a function of the velocity at the 
NLO, when the thermal mass of a parton is 400 MeV(left) 
or 600 MeV(right) respectively. These are the sum of 
quark induced and gluon induced NLO processes. As 
in the LO case, the width is smaller for larger thermal 
parton mass. This is mainly due to a larger virtuality 
of the gluon propagator when the initial parton has a 
larger thermal mass. As shown in the figure, in contrast 
to the results at the LO, the thermal widths slowly in- 
crease as the velocity of J/tp increases. Such different 
behavior results from the asymptotic form of the NLO 
dissociation cross section, which converges to some finite 
value at the high energy limit. In the appendix, we de- 
rive the asymptotic form of the NLO dissociation cross 
section, assuming the thermal parton mass to be suffi- 
ciently small. If the cross section is almost constant, the 
thermal width can be approximated as follows, 



T - a 



d 3 k 

(27T) 5 



n(k )v re i((3) ~ a < v rei (f3) >, (5) 



where < v re i(j3) > is the expectation value of the rela- 
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FIG. 2: The variation of the thermal width T of T in the 
LO (upper) and in the NLO (lower) as a function of their ve- 
locity with the parton thermal mass of 400 MeV(left), of 600 
MeV(right) 



process. This causes the monotonically increasing NLO 
cross section to reach its asymptotic value at a higher 
incoming energy when the binding energy is larger. This 
is the reason why for the T system, the thermal width 
has a larger slope when the velocity increases. The ther- 
mal width steadily increases until the T is traveling close 
to the speed of light so that most of the thermal par- 
tons have sufficient amount of energy to dissociate the 
T at the asymptotic limit. In contrast, for the J/tp at 
rest, non trivial fraction of the thermal partons already 
dissociates the quarkonium at the asymptotic limit. 

Up to now, the thermal mass of partons were assumed 
to be temperature independent for simplicity. However it 
is predicted to scale as g(T)T by finite temperature QCD 
calculations. Therefore, we finally present the result for 
the thermal widths of J '/ip and T obtained with tem- 
perature dependent thermal masses for the partons. The 
masses of thermal gluon and quark are taken respectively 
as, 



with 



m 2 g (T) 
m 2 q (T) 



9 2 (T) = 



g 2 {T)T 2 



9 2 {T)T 2 



3 



N F 

IT 



(G) 



48tt 2 



(HiV c -2iV / )lnF 2 (r,T c ,A)' 



The number of color N c is set to 3, and the number of 
flavor Nf to 3. For the function F(T,T C , A), we use the 
form obtained in from a fit to the lattice QCD calcu- 
lations. Fig. [3] shows the thermal widths of J/tp and of 
T obtained with temperature dependent parton masses. 



tive velocity between partons and the J/tp moving with 
velocity (3. As mentioned above, as J/tp moves faster, 
the expectation value of the relative velocity increases. 
If the J/tp moves with the speed of light, < v re i((3) > 
approaches 1, and r ~ a. In conclusion, while the width 
decreases in the LO, its magnitude is small such that the 
sum of the LO and the NLO width increases with the 
velocity. 

Fig. [5] shows the thermal widths of T as a function of 
its velocity in QGP. As in the case of J/tp, the widths 
decrease in the LO, but increase in the NLO as the ve- 
locity increases. However, the slope of increase is steeper 
for T than for the J/tp, which comes from the differ- 
ence in their binding energies. If the binding energy is 
larger, the threshold energy becomes larger and so does 
the energy at which the LO cross section becomes maxi- 
mum. The dominant NLO process is a forward scattering 
contribution, where the incoming parton emits a virtual 
gluon, which then dissociated the quarkonium via the LO 
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FIG. 3: The variation of the thermal width F of J/tp (left) 
and of T (right) from the sum of LO and NLO as a function 
of their velocity, obtained with the temperature dependent 
thermal mass [401 ] 
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The magnitudes of the widths lie between the boundaries 
obtained with a constant thermal mass of 400 MeV and 
600 MeV, as shown in Fig. Q]and in Fig. [2]for J/ip and T 
respectively. This is so because the temperature depen- 
dent thermal mass lie between 400 MeV and 600 MeV 
in the temperature range considered in this work. More- 
over, the functional dependencies of the thermal widths 
with respect to the velocity are similar to those obtained 
with a constant thermal mass. 



IV. DISCUSSION 

In the heavy ion collisions at LHC, not only will the 
heavy quarkonia be more amply produced, but be pro- 
duced more energetically. The suppression and/or rel- 
ative enhancement of these high pr quarkonia will also 
occur in such an environment. In this respect, it is very 
important to know how the properties of heavy quarko- 
nia will change when they move with respect to the QGP. 
Our result, based on the NLO perturbative QCD calcu- 
lation, shows that the thermal width becomes larger as 
the quarkonium travels faster with respect to the QGP. 
The rate of increase is larger for the case of T case than 
for the J/ip, although the magnitude itself is larger for 
the latter. Our result suggests that the survival rate of 
J /if) will be lower than that obtained with the velocity 
independent thermal width calculated at restJUJ. Con- 
sequently, if a relative suppression of J/ip with higher 
than lower pt is observed, as suggested by (26|, it may 
be a consequence of the broadening of the thermal width 
of J/tp as well as of the shortening of the Debye screen- 
ing length. The difference between these two effects is 
that the former depends on the size of fireball while the 
latter does not. Therefore, a systematic study on the 
A-dependence will be able to discriminate between these 
two effects. 



where q, k±, ki are the momentum of quarkonium $, 
incoming thermal quark, and outgoing thermal quark re- 



spectively, and p\ — (fci — k 2 ) 



w 



(q + PA) 2 - The 



momentum p is the relative three momentum between Q 
and Q. If the quarkonium is the Coulomb bound state 
of 15, the absolute square of derivative of quarkonium 
wavefunction is 



dip(p) 



dp 



«oP 2 
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2 10 7r(a e 



1)6 

- k 2 o 



(8) 



where ao is the Bohr radius and eo is the binding en- 
ergy of the quarkonium. In the second line of the above 
equation, the energy conservation condition to$ + ki$ — 
P 2 / TO Q + ^20 j ancl the relation a 2 , = l/(e mg) have been 
used. It might seem odd to use the nonrelativistic energy 
conservation condition in the high energy limit. However, 
in the high energy limit, the forward scattering is domi- 
nant, and therefore, while the incoming quark energy is 
very large, the transferred energy carried by the virtual 
gluon that dominantly dissociates the quarkonium is not 
large. This is so because the LO dissociation cross sec- 
tion is dominant near threshold. Hence a nonrelativistic 
treatment of the bound state kinematics is justified, a, 
/?, a', and (3' are the integration limits for drawing Dalitz 
plot on the plane of p\ and w . They are respectively 



a = 4m; 
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= (y/s — m kl ) 2 
a' = —b— \/b 2 — ac 
P = -b + \Jb 2 - ac, 



(9) 



where 
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Appendix 



{s - (m* + m kl ) 2 }{s - (rn$ - m kl ) 2 }/(2s) 
-{s - (ml - m 2 kl )}(w 2 - m%)/(2s) 
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x{w 2 - (u + m kl ) 2 }{w 2 - ( M -m fcl ) 2 }/(4s 2 ). 

Here, s is the square of the initial energy in center-of- 
mass frame. In the limit of s — > oo, 



Here, we derive the asymptotic form of the cross sec- 
tion for the process & + q^Q + Q + q in high energy 
limit. The cross section for three-body decay is expressed 
as follow 1311 
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and the initial flux becomes 
4 



l^J(q ■ kx) 2 - m 



2 9 



2s. 



(10) 



6 



Then the elementary cross section becomes 



16g 4 TOgTO$(aoeo) 
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Because k\$ — A^o varies from eo approximately to s 1 / 2 , 
first two terms in the square bracket have no contribution 
in the large s limit. After integration with respect to p\, 
the cross section becomes 



As the forward scattering becomes dominant in high 
energy scattering, most of the contributions come from 
the region w 2 ~ 4toq in the integration of w 2 . There- 
fore, the the main suppression factor in the integrand in 
Eq. (fl"2|) in the s — > 00 limit is l/(w 2 — m|) 5 . Terms 
such as ln(u; 2 — m|,) and l/Vw 2 can be replaced by 
ln(4mg — m\) and 2toq respectively. Then we have, 
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where we used the following decomposition formula for 
integration with respect to p\ , 
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(a — b) n x — b (a — b) n x — a 

and we ignored the thermal mass for simplicity ex- 
cept in the logarithm as a regulator. This elimination 
brings about a potential problem when the binding en- 
ergy is very small. To see its subtlety, consider the de- 
composition formula in the following example, 
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Within the integration range, the left side of Eq. (TT3"]) is 
always positive. However, if the binding energy is very 
small, l/(w 2 — m| — 2m\ ) of the right side of Eq. (TiU|) 
has a singularity. Of course, the term in the parenthe- 
sis on the right side diverges at the same time, and the 
equality is maintained. However, if the thermal mass 
2m 2 i is dropped during the approximation for simplic- 
ity, such divergence might not cancel and bring out the 
wrong result. 
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the w 2 dependent terms become 
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and finally the cross section becomes 
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Fig. |4] shows the discrepancy between the result of Eq. 
and the result of Eq. (T5j) at y /s = 100 GeV, where 
we used the binding energy of J/i/j in the vacuum [35j ] - 
It is shown that the error decreases as the thermal mass 
of a parton decreases. If we further ignore the binding 
energy of a quarkonium, m$ = 2mQ, we find 



48tt 



1)17 



97 
20 



(16) 



This formula has several important features. First, the 
cross section is proportional to the square of the Bhor 
radius. It is a consequence of the dipole type of cross 
section, which is proportional to the derivative of the mo- 
mentum wave function squared. The second important 
feature is the logarithmic term with argument propor- 
tional to the thermal mass, which acts as the regulator. 
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FIG. 4: The discrepancy between Eq. (J7J) and Eq. (fTo)) 



In the case of $ + 5— * Q + Q + g, the terms in the 
square bracket of Eq. (fTTj) are replaced by the following 



2k 



10 



2 A: 



20 



V20 



>10 



Z-2 

"-20 

^10 



k 2 

k 2 



fci • fc 2 



(fcio - fc 2 o) 



2 J (^10 



- 2 



(17) 



In the asymptotic limit, all but the last term in the 
square bracket of Eq. (fT7|) are suppressed in the large 
s limit. Because its coefficient is twice of that for the 
quark induced case, the asymptotic cross section is two 
times larger. This ratio can also be seen in the quark and 
gluon induced dissociation cross sections shown in figures 
of ref [3l|. 

We remark that the asymptotic value for the cross sec- 
tion obtained in the Appendix is for the coulomb bound 
state, while in the previous sections, both the Bohr ra- 
dius ao and the binding energy eo of quarkonia are ex- 
tracted independently from lattice data, such that the 
relation a 2 , = 1/(£otoq) is not satisfied. Moreover, the 
extracted binding energy is too small compared to the 
thermal mass of a parton, especially in the case of J /ip. 
Because of these reasons, Eq. flTSj) or Eq. ([To]) cannot be 
used directly. Nevertheless, the property that the disso- 
ciation cross section in the NLO converges to nontrivial 
finite value in the high energy limit is still maintained. 
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